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Renormalization Group Approach to Oscillator Synchronization 
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We develop a renormalization group method to investigate synchronization clusters in a one- 
dimensional chain of nearest-neighbor coupled phase oscillators. The method is best suited for chains 
with strong disorder in the intrinsic frequencies and coupling strengths. The results are compared 
with numerical simulations of the chain dynamics and good agreement in several characteristics is 
found. We apply the renormalization group and simulations to Lorentzian distributions of intrinsic 
frequencies and couplings and investigate the statistics of the resultant cluster sizes and frequencies, 
as well as the dependence of the characteristic cluster length upon parameters of these Lorentzian 
distributions. 

PACS numbers: 



I. INTRODUCTION 

One of the most exciting aspects of modern physics 
is the investigation of emergent collective structures in 
many-body problems. Spontaneous synchronization in 
networks of interacting nonlinear oscillators with differ- 
ent (often randomly distributed) frequencies is one of the 
best understood collective phenomena in noncquilibrium 
systems If the interaction between the oscillators 
is weak, then the oscillators will evolve in an uncoordi- 
nated fashion. In contrast, with proper connectivity and 
coupling strength, the population can display the same 
frequency as order forms. This collective behavior has 
been observed in diverse systems including neural net- 
works @, H, 0|) Josephson junctions @, @, 0|, lasers [|[, 
electronic circuits @, and a range of biological systems, 
including sleep cycles [loll, animal gaits [1JJ and biologi- 
cal rhythms in general [12j . 

Useful insights into the phenomenon of synchroniza- 
tion are obtained by analyzing the idealized case of all- 
to-all coupling in the limit of a large number of oscil- 
lators. In one of the most common models each oscil- 
lator is described by a single variable, phase, that ad- 
vances at a constant rate in the absence of coupling. This 
model is exactly soluble using a mean-field approach, and 
it provides a framework for understanding synchorinza- 
tion [ID, 0, [H[ • In this model, above a critical coupling 
strength the synchronized phase appears: a finite fraction 
of oscillators evolve at a common frequency. Many fea- 
tures of this phase, such as the growth of the fraction of 
locked oscillators as a function of the coupling strength, 
can be calculated. 

In many physical systems finite range interactions are 
a more realistic description. The simplest of these inter- 
actions are nearest neighbor couplings [HI, [TlJ • In gen- 
eral, finite range coupled systems are much harder to 
analyze, and, as a result, there are few theoretical tools 
that produce a quantitative description of the collective 
motion. Much of our understanding of these systems 



rests on numerical solutions of the governing equations. 
We propose an alternative method to attack the synchro- 
nization problem. The method we develop is a real-space 
renormalization group (RG) designed to work well on sys- 
tems with strong randomness, where both the intrinsic 
oscillation frequencies and coupling constants are ran- 
dom variables taken from distributions with long tails. 
The specific example we investigate contains Lorentzian 
distributions. This method is shown to be accurate in 
systems with short range interactions. 

In this paper, we analyze a one dimensional chain of 
oscillators interacting with their nearest neighbors. De- 
spite the relative simplicity of the model, it can exhibit 
complex behavior since randomness competes with a ten- 
dency to establish macroscopic order. This system is 
below the lower critical dimension for a synchronization 
transition (Strogatz and Mirollo [l8|, [l9| showed that a 
nearest-neighbor chain cannot exhibit extensive synchro- 
nization), and therefore can not be analyzed with mean- 
field techniques. Collective behavior still arises, albeit at 
a finite length scale: the chain becomes fragmented into 
finite length clusters of coherently moving oscillators. We 
are interested in the distribution of cluster lengths, clus- 
ter frequencies, and how these change as a function of 
the coupling strengths and bare frequencies. Addition- 
ally, we seek to develop a renormalization group scheme 
that will allow us to predict with some accuracy the in- 
dividual local behavior of oscillators. We also carry out 
an extensive numerical study of the problem and use it 
as a benchmark for comparison with the RG. 

The paper is organized as follows. In Sec. [TT]we outline 
the strong disorder renormalization group approach with 
the details of the steps for the nearest-neighbor oscilla- 
tor chain derived in Sec. Ill Al In Sec. IIII1 we compare 
some key statistical features of the cluster state obtained 
by the RG with results from direct numerical simula- 
tions of the equations of motion. In Sec. IIVI we outline 
our findings concerning the physics of the synchronized 
clusters, presenting results for the distributions of cluster 
sizes, cluster frequencies, as well as the dependence of the 
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characteristic length scale of the clusters on the widths 
of the coupling and frequency distributions. Concluding 
remarks are made in the last section. 

II. STRONG DISORDER RENORMALIZATION 
GROUP 

To analyze the random oscillator chain we develop 
a strong-disorder real space RG technique. A similar 
method has been used successfully in quantum mechan- 
ical systems to analyze the ground state of the random 
Heisenberg model [2(| [U [22|, HH, and the superfluid- 
insulator transition in random bosonic chains [13, HH (a 
quantum reactive analog of the nearest-neighbor chain, 
where the 9 replaces 9, and 9 is a quantum operator). 

The main idea behind the strong disorder approach ap- 
plied to the quantum-mechanical problem is best demon- 
strated in the case of the spin-1/2 Heisenberg model: 

H = J2 J iSi-Si+i- (!) 

i 

The i'th term in the sum considered on its own, splits the 
four states of the spins i and i + 1, into a singlet state, and 
three triplet states excited by the energy J^. Thus a real- 
space decimation step can be proposed. First, freeze the 
spin pair with the strongest interaction, J n = max{Ji}, 
into a singlet state. Perturbative corrections to this state 
introduce an effective coupling between the neighbors of 
the spins n and n + 1 which is: 

tt _ Jn-lJn+1 g ,2, 

tlRG — ^TJ Es n -a n +i, (2) 

an interaction identical in form to the operators in the 
bare Hamiltonian, but with a much suppressed strength: 

Jeff = J "~2jl" +1 < Jn-i, J n , J-n+i- By repeating this 
decimation step, the number of free spins is gradually 
decreased; at the stage when all spins are bound into 
singlets, we obtain the ground state of the Hamiltonian. 

The decimation step described above is justified if the 
coupling to the rest of the chain is much weaker then 
between the two spins in question. This is guaranteed if 
the chain considered has broad coupling distributions. 
Below, we use a real-space renormalization group ap- 
proach to analyze the nearest-neighbor oscillator chain 
with broad disorder distributions. 

The nearest-neighbor oscillator chain presents a classi- 
cal problem with dissipative equations of motion, rather 
than a quantum problem with a conserved Hamiltonian. 
Nevertheless, we can develop an analogous strong disor- 
der RG method for this problem as well. The chain is 
governed by the equations of motion 

6i = Ui + K t -x sin (6>i_i - 9i) + K t sin (9 i+1 - 6i), (3) 

where 0j is the phase of the ith oscillator. The uji are 
the intrinsic frequencies taken from a random distribu- 
tion (we assume zero mean without loss of generality) and 



the Ki give the couplings to the nearest neighbors also 
taken from a random distribution. The coupling will or- 
ganize the oscillators into clusters of common frequency 
Oj defined as 

^ km ftW-W (4) 

(t— to)— >oo t — £q 

such that each oscillator in a given cluster shares its Q 
with all the other oscillators in the same cluster. We 
seek to understand the statistics of the cluster sizes and 
frequencies. 

The renormalization group approach is based on two 
observations. A strong coupling between two oscillators 
tends to force them to rotate at the same frequency form- 
ing a cluster that we can then treat as a single effective 
oscillator. Concomitantly, oscillators with anomalously 
large frequencies rotate essentially independently, with 
little effect from their neighbors. These ideas lead to two 
decimation steps: a strong coupling decimation and a fast 
oscillator decimation. We apply these steps iteratively, 
sweeping through the oscillators from large values of u>i 
and Ki to smaller values. 

A. Decimation steps 

As in the case of the quantum Heisenberg chain, our 
real space RG will be carried out iteratively by treating 
at each step the strongest term in the right-hand-side of 
the equations of motion, i.e., the term with the highest 
frequency. In the oscillator chain the strongest term at 
each step could be either a strong coupling Ki, or a large 
frequency, Wj. Let us now derive the appropriate deci- 
mation steps for the nearest-neighbor chain for these two 
cases. We begin by slightly generalizing the model to 

rriiOi = m^i + Ki-i sin - 9 { ) + K, t sin (9 i+1 - 9, t ), 

where a new parameter m is added. It will follow from 
the strong coupling decimation step defined below that 
when two oscillators with rrii and m^+i are described by 
one effective oscillator, its m is given simply by mi+mj+i. 
Hence, the value of m represents the number of original 
oscillators that this effective oscillator represents. It will 
become apparent that oscillators with larger m are less 
affected by perturbations from their neighbors, so m is a 
measure of the inertia of an oscillator and will be referred 
to as its "mass" , although it is not a true mass or inertia 
in the sense of Newton's second law. 



1. Strong coupling decimation 

The phases of two oscillators connected by a very large 
coupling tend to advance in synchrony. The strong cou- 
pling step replaces this pair of oscillators by a single ef- 
fective oscillator with a different mass and frequency de- 
termined as follows. Consider the largest coupling in the 
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FIG. 1: A small piece of the chain around the strongest cou- 
pling, labeled here as K n . All the neighboring Ki and all the 
local u>i are assumed to have much smaller magnitude then 



chain: suppose this coupling is K n , Fig.[TJ Because of the 
strong disorder the frequencies oj n , w n +i of the pair of os- 
cillators that K n couples and the couplings K n -i, K n +i 
to the neighbors of the pair will almost always be small 
compared to K n . Ratios such as oj n /K n , K n -i/K n , are 
therefore expected to be small numbers; we will use a 
symbol e to denote quantities of this order. For small e, 
K n will couple oscillators n and n + 1 into a synchro- 
nized cluster that can be represented as a single effective 
oscillator described by phase 

6 = (m n 6 n + m n+1 e n+1 )/M, (6) 
possessing mass M, 

M = m n +m n+1 , (7) 
and an effective intrinsic frequency f2, 

SI = {m n uj n + m n+ iL0 n+ i)/M. (8) 

We refer to such pair of oscillators strongly coupled 
(cf. H| and Sec. HEB]). 

The small frequency difference u) n +i —u> n and t ne small 
couplings K n -i and K n+ i induce a small phase difference 
S = n +i — n between the two oscillators. Using 9 n — 
— m n +\8/M and 9 n+ i = + m n S/M, Eqs. (O gives 
the evolution equation for the phase of the new effective 
oscillator 



A/0 = Mil + sin (0 n _ a - 



M 



5) 



+ K n+1 sm(9 n+2 -e-j^5). (9) 

There are corresponding equations for 9 n —i and 6 n +2- 
Aside from the new <5-terms, the equations take the orig- 
inal form. 

There are two contributions to the phase 8. The con- 
tribution from the frequency difference is a time indepen- 
dent phase given to leading order in e by 



/^n,n+l 



n+1 



' W„). 



(10) 



where ^ = m- 1 



_1 is the reduced mass. This can 



be seen by considering just the two oscillators coupled by 
the strong K n and disregarding effects of the rest of the 
chain. Then, as can be seen from the difference of their 
equations of motion, Eqs. ([5]), as long as w n +i — u> n is 



smaller then K n J fi n>n +x, which will be the case if K n is 
much larger then all the neighboring parameters, these 
two oscillators will reach a state in which they rotate at 
the same frequency and are separated by the fixed phase 
difference <5o- We will see this explicitly after writing 
down an equation of motion for <5 below. For the one di- 
mensional chain we can redefine 6 n —i and prior phases, as 
well as 9 n+ 2 and subsequent phases, to absorb such con- 
stant phases, and so we may ignore them. This maneuver 
is restricted to one dimension in the spatial lattice, and 
the extension of the method to higher dimensions would 
require a more careful study of the effect of these addi- 
tional phases. 

The couplings K n -\ and K n+ i to the neighbors of the 
pair give an additional contribution to 6 that is time de- 
pendent as the neighboring oscillators evolve. To leading 
order in e it turns out that these oscillating corrections 
give small extra coupling terms of forms not included 
in the equation of motion Eqs. a second nearest 

neighbor interaction between oscillators, a three body 
interaction involving the new effective oscillator and the 
neighboring pair, and a second harmonic correction to 
the form of the pairwise interaction. The detailed forms 
of these terms can be found in [2(| . We will neglect these 
more complicated interaction terms since they are much 
weaker then all other competing interactions. 

These results arc derived from the exact expression for 
the equation of motion of the phase difference S(t), 



S = — 



Kn 



Mn,n+1 
K n +1 

m n +i 



■ Sin (5 + {ljJn+1 - CJn) 



• sm 



'n+2 







■ sm 



h-i - 



M 
m n+ i 
M 



(11) 



This equation describes a damped particle in a time 
dependent washboard potential. The relaxation rate 
^n//*n,n+i given by the first term on the right hand side 
is fast compared with the other terms in the equation. 
The solution can therefore be developed as an expansion 
in a ratio of time scales, and the leading order term is 
the adiabatic solution 



fJ"n,n+l i \ . ^n+1 . , 

— - — (w„+i - u n ) H sm ( 

Kn Y m n+ i 



'n+2 



-0) 



'n-\ 



0) 



(12) 



The first term in the braces gives the constant phase Sq, 
and on substituting back into Eq. ^ and expanding to 
first order in small (5, the second and third terms give the 
corrections to the interaction terms that we will neglect, 
as discussed above. Specifics of the calculations outlined 
in this section are presented in [26| . 

In summary, the effect of the strong coupling decima- 
tion step is to replace the pair of oscillators with phases 
6 n and 9 n +i by a single oscillator with phase 0, mass M 
and frequency O given by Eqs. ©, and ©. 
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FIG. 2: A small piece of the chain around the oscillator with 
the largest frequency, labeled here as uo n . All the neighbor- 
ing Ki and all the neighboring Ui are assumed to have much 
smaller magnitude then u> n . 



2. Fast oscillator decimation 

Next we consider the case where the highest frequency 
scale in the problem is a bare frequency. In this case 
we apply the notion that the phase of an oscillator with 
frequency much different from its neighbors tends to ad- 
vance as if the oscillator were decoupled. Consider the 
oscillator with the largest frequency in the entire chain: 
suppose this is oscillator n with intrinsic frequency uj n , 
Fig. O Because of strong disorder the coupling strengths 
K n -i and K n to the neighboring oscillators and the fre- 
quencies uj n ±i of these oscillators will almost always be 
small compared to u> n . Therefore we expect ratios such 
as K n _i/uj n and u) n -i/u) n to be small compared to 1; as 
before we denote quantities of this order by e. To ze- 
roth order in e the n'th oscillator runs freely, with the 
frequency ui n . We call such oscillators "fast oscillators" 
(see and Sec. Ill Bl for elaboration). 

Deviations from the free-running solution can be cal- 
culated perturbatively in e. The most important effect 
we might look for is an induced effective coupling via os- 
cillator n between the neighbors n — 1 and n+1, since this 
might lead to their synchronization. However, at least to 
order e 2 as we describe below, we find that there is no 
induced interaction that tends to lock the neighboring 
oscillators 31]. Thus the main effect of decimating the 
largest frequency oscillator in the chain is to eliminate 
this oscillator from further consideration (it represents 
a cluster in the final tally at frequency close to uj n and 
with a size given by its mass parameter m) and to cut 
the chain into two pieces. This is a key feature of the 
strongly random chain that limits the size of synchro- 
nized clusters. At order e 2 we do find a renormalization 
of the frequencies of both the eliminated oscillator and 
the two neighbors remaining in the chain. 



We now investigate the perturbation of the motion of 
the fast oscillator (oscillator n in Fig. [J) given by Eqs. (JSJ). 
To zeroth order in e, the solution to n is given by the 
uniform running state n — u) n t. There are two 0(e 2 ) 
effects that slightly perturb the average rate of this uni- 
form advancement of the phase. The first is the fact that 
oscillator n rotates in the tilted washboard potential cre- 
ated by its neighbors. Since typically o;„_i,w n +i -C w n 
for strong disorder, this potential can be treated as ef- 
fectively static. The second effect is a small component 
of fast dynamics added to the motion of the neighboring 
oscillators, which then acts back on oscillator n. Both 
of these effects change the average frequency of the elim- 
inated oscillator by a shift of order e 2 . There are cor- 
responding corrections to the motion of the neighboring 
oscillators as well. These effects will be apparent in the 
following analysis. 

First, consider the motion of n in the nearly static 
washboard potential from the neighbors, calculated to 
second order in the coupling. Setting the neighboring 
phases in Eqs. (JSJ) to fixed values 0^±i gives the oscilla- 
tion period of the nth oscillator. 



T„ 



d6r. 



o uj n + sin 10 



Expanding the denominator and performing the integra- 
tion over 9 n gives the change to the frequency of the nth 
oscillator Q n — ui n with ui n — 2irT~ 1 



K 2 



K n -\K n 



K 2 

2m 2 w„ 2m 2 l u n m 2 n uj n 



q(0) 

y n+l 



j(0) 



(14) 



The second effect perturbing the motion of the fast 
oscillator arises from the rapid induced motion of the 
neighboring oscillators n ±i — t } l ±i+S0 n ±i- For example 
from Eqs. © for oscillator n — 1 we have 



50n- 



K n - X 
rrin-i 



m„_if2 



dt sin (uj n t — 0^-1) 



i-COS (w n t -<?£_!), 



(15) 

(16) 



where we use the zeroth order solution for oscillator n 
and the influence of the next nearest neighbor 9 n -2 only 
appears at higher order. This provides a perturbation 
back on the nth oscillator given by expanding the second 
term on the right hand side of Eqs. [5] for i = n to first 
order in 50 n -i 



S0 r , 



Kn-l 



COS (u> n t 



~>n-l- 



(17) 



Using Eq. (fT5|) for 50 n -i, averaging over the fast oscilla- 
tions, and adding the similar effect from oscillator n+1 
gives the change in the average frequency of the nth os- 
cillator 



K 2 



K 2 



2m„m„_iw„ 2m n m n+ iuj n 



(18) 
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Combining Eqs. (fT4|) and (fT8|) gives the total renormal- 
ization of the frequency of the fast oscillator to second 
order 



SLOn 



K 2 



K 2 



K n -\K n 



m 2 uj n 



■ cos ( 



(19) 



m ■ 1 as before. 



with ^r. = m" 1 

There are reciprocal effects on the neighboring oscilla- 
tors that can be calculated by a similar procedure. Since 
at the pairwise level interactions between the oscillators 
conserve the mass- weighted frequency average [32| . we 
can immediately see that the first two terms in Eq. (fTO)) 
correspond to frequency renormalizations of the neigh- 
boring oscillators 



Suj n - 



K 2 



2m r , 



-lfJ-r. 



K 2 



+ 1 



(20) 

The fact that the first two frequency correction terms 
in Eq. (|19p are additive suggests that they can be ob- 
tained from the analysis of a pair of 2-oscillator systems. 
This calculation can be performed exactly for the phase 
difference, for example <f> = 9 n — 6 n -i- Such a frame 
independent analysis reproduces at 0(e 2 ) the first two 
terms in Eq. (fl9|) and the terms in Eq. ([20|) which were 
obtained in the specific reference frame of the average 
frequency, but also yields terms that are higher order in 
e. In the context of the full chain including these higher 
order terms is an uncontrolled approximation, however 
empirically we find them to produce a better match with 
the simulation data in some rare cases. Hence, in the 
numerical renormalization of the chain, we employ the 
expressions produced by this analysis. 

The phase difference of oscillators n and n — 1 isolated 
from the rest of the chain satisfies 



(u n - uj n -i ) sin ( 

^n,n— 1 



(21) 



A period T over which this phase difference grows by 27r 
can be defined analogously to Eq. (JT3J), but since there 
is only one sine term, the integral can be computed an- 
alytically. From that solution and the fact that for the 
isolated pair m\Q\ +m.2#2 = m\(j)\ +rri2Uj n we obtain the 
correction to the frequency of oscillator n coming from 
oscillator n — 1: 



Mn,n-l(^n-l 



Kn-1 



Mn,n-l(k>n — U> n -i) 



The new cosine interaction term in Eq. (| 19[) also war- 
rants discussion. There is a corresponding renormaliza- 
tion of the equations of motion of each of its neighbors. 
In the motion of the fast oscillator Eq. (JT3J) , we ignore 
the cosine term since it averages to zero on the long time 
scale of the period of oscillators nil. In the equations 
for # n _iand 9 n+ \ the corresponding term is a coupling 
oc cos(#„ + i — (9„-i) between these oscillators, with the 
same sign in the equations of motion of both neighbor- 
ing oscillators. Therefore this interaction does not have 
action-reaction symmetry. This sign peculiarity arises 
because the decimated oscillator has a particular direc- 
tion of rotation. Because of the lack of action-reaction 
symmetry, the cosine term does not tend to pull neigh- 
boring oscillators together, and so does not directly in- 
fluence the clustering. We therefore omit this term in the 
renormalization procedure. 



B. Implementation 

The chain of oscillators is renormalized by successive 
application of the two decimation steps developed above. 
These are executed numerically on a list of parameters 
(mj, u>i, Ki) representing the chain of oscillators. A single 
decimation step consists of finding the largest frequency 

of the set iu>i : 2 ^ <i +1 y renorma li z i n g the appropriate 
term as described above, and working down to smaller 
ones later. Notice that it is the largest tt-^ 1 — that iden- 
tifics the pair of oscillators to be subjected to the strong 
coupling step, not the largest Ki. The reason for this 
becomes clear when all equations are divided by their re- 
spective mi . An equation for the phase difference 0j+i— 6 
contains the term 



sin (0j+i —9i). When 



Ki 



is 



(22) 



We use this expression to replace the first term on the 
right hand side of Eq. (|19p with a corresponding expres- 
sion for the first term in Eq. (|20[) . and similar expressions 
for the second pair of terms. 



Hi,i+i ~ v " '~ ri " •'" ■ ■ /*<,«+ 1 

much larger then the terms representing coupling to the 

neighbors, this phase difference will tend to be locked and 

the strong coupling step should be applied to the pair of 

oscillators i and i + 1. The factor of 2 was used in the 

denominator to make sure that the initial stages of the 

RG procedure, couplings and frequencies are treated on 

equal footing. 

Initial values of masses are all set to 1 . The numerical 
procedure identifies the oscillators with to and K that 
lie in a narrow band of magnitudes at the top of the 
spectrum of remaining values, and then decimates those 
oscillators according to the steps defined in Sees III A II 
and III A 21 A band is used, rather than just selecting 
the largest values, to improve the efficiency of the code. 
The width of the band is chosen to be narrow (1% of 
the chain) to maintain the descending order in lu and K 
for the decimation of nearby oscillators. Oscillators dec- 
imated in the fast oscillator step form a cluster in the 
final tally, with the effective mass describing the size of 
the cluster and the frequency giving the time-averaged 
rate at which the phase of this cluster advances. As de- 
scribed in Sec. Ill A 21 these oscillators are removed from 
the chain. After the decimation of oscillator band, the re- 
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mainder of the chain consists of fewer oscillators, which 
have lower coupling strengths and frequencies within a 
narrowed spectrum. The procedure is then repeated for 
successively lower bands of K and u>. 

During the RG execution, two types of special cases 
related to violations of the assumptions based on strong 
disorder are occasionally encountered. First, in our treat- 
ment of the oscillator with the largest u> we have made an 
assumption that nearby u> and K are small in compari- 
son. This assumption defines the physical content of the 
strong disorder assumption. In real chains, it is possible 
to have neighbors with parameters which are somewhat 
lower yet the difference is not large enough to guarantee 
that effects of these neighbors are small. This is espe- 
cially important as the spectrum of the remaining K and 
u> values shrinks. To be more precise, for a decimation of 
a fast oscillator n, the magnitude of rpo defined by 



tfo = Max 



Kn 



(23) 

determines whether the effects of the neighbors are small, 
so that the step can be performed. Similarly, for a strong 
coupling decimation of oscillators n and n + 1, the mag- 
nitude of the ratio 



K„ 



rsc 



Hn,n+l\<^n — ^n+l| 



(24) 



determines whether the step should be performed. In 
our implementation of the RG we set set the criteria to 
rpo < 1 and rsc > 1- This choice was based on the 
numerical study of small chains (N < 10) as well as on 
the fact that for sinusoidal coupling t F q does not need 
to be very small for uj n t to give a good approximation 
to the dynamics of the nth oscillator. This is hinted by 
the square-root dependence of the frequency correction 
in Eq. (f2"2"|) : analysis of that formula indicates that ui n 
quickly approaches u) n as Tfo decreases below 1. A more 
thorough discussion of the choice for rpo an d r sc can 
be found in [26| . 

There is a second type of violation of the assumptions 
based on the strong disorder. Note that single oscilla- 
tors with Tfi — 1 or oscillators representing clusters with 
to > 1 are normally subjected to the fast oscillator deci- 
mation step at some point in the RG procedure: most of 
the bare oscillators either join other oscillators in build- 
ing clusters, which are eventually decimated via the fast 
oscillator step, or they are decimated as fast oscillators 
directly without undergoing a strong coupling step. But 
some rare combinations of frequencies and couplings will 
prevent the decimation from occurring at all within this 
scheme. For example, while an oscillator's frequency may 
lie in the executable band, this oscillator may be coupled 
via a K which is outside this band (it has a somewhat 
lower magnitude), but may nevertheless cause rpo to ex- 
ceed 1. In general, when the RG reaches the scale of that 
coupling, it will perform a strong coupling decimation 
step. However, in rare cases this may not happen due 



to intervening steps that have modified the neighboring 
parameters, rendering that coupling no longer strong ac- 
cording to the rsc criterion. As a result, the site in ques- 
tion has not been subjected to either decimation step. To 
address such unusual cases, the algorithm sweeps repeat- 
edly through the entire spectrum of ui and K values until 
nothing is left to decimate. 



III. NUMERICS 

While the discussion above constructed a general RG 
algorithm, let us now explore the scheme within a dis- 
order realization with a specific class of frequency and 
coupling distributions g(w), G{K). The bare frequencies 
in the system we analyze have a symmetric Lorentzian 
distribution with a cutoff: 



i 



(25) 



with —uj c < u> < lo c . The couplings are positive and are 
taken from a half Lorentzian with a cutoff: 



G{K) 



2/jC 2 



K 2 



(26) 



with < K < K c . The cutoffs are used to facilitate 
the effective numerical simulations fSec. lIIIBl below): in- 
creasing cut-off values increases the time needed for sim- 
ulations. The same cut-off values for {w} and {K} were 
chosen: us c = K c — 100; the values of these cutoffs arc 
much larger then the FWHM of all Lorentzians that we 
considered. The C\ and C 2 are normalization constants 
that become 1 /it when the cut-off values uj c and K c are 
infinite. Notice that due to the structure of Eqs. (O, it is 
always possible to divide all equations by a constant such 
that the width of one of these distributions will be unity: 
this only changes the time scale of all the processes in the 
chain. We chose g(u>) to have unity width and explore 
the physics by varying the other width parameter /i. 



A. Renormalization group 

The RG procedure is applied to a chain of 10 6 oscil- 
lators. Twelve values of coupling width /i are studied: 
jU = 0.025, 0.05, 0.25, 0.5, 0.625, 1.25, 2.5, 3.75, 4.5, 5, 
6.25, and 7.5. The set of random numbers defining the 
particular realizations of the chain were different in the 
RG and in the simulations (below) for the statistical com- 
parisons, but were identical for the real-space comparison 
described in Sec. IIII CI 



B. Simulations 

In order to explore both the system's characteristics 
and the reliability of the renormalization group, Eqs. (J3J) 
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are also integrated numerically. The numerical method 
is a variable stepping Runge-Kutta algorithm. Systems 
of N = 10000 oscillators are solved with N intrinsic fre- 
quencies chosen randomly from the cutoff distribution 
given by Eq. (|25|) . Similarly N — 1 coupling constants 
are randomly selected from the cutoff distribution given 
by Eq. (126|) . The same twelve values of coupling width fi 
that are used in the RG are studied in the simulations. 
Results for each value of \i are averaged over 100 realiza- 
tions of the intrinsic frequencies and coupling constants. 
To help facilitate comparisons between different distribu- 
tion widths the same 100 random number seeds are used 
for each distinct coupling width. 

The simulations are performed for a relatively long 
time of t — 10000. Running phases at each site are 
recorded at regular intervals and used to calculate av- 
erage frequencies over a time T according to 

*<*> - fl '" 0+ V"° > ; < 27 > 

with T = t — <o [compare with the theoretical definition 
of Q given by Eq. The time to is chosen to elimi- 

nate transients: increasing values are tried until average 
frequencies remain unchanged. In these simulations, the 
value of 2tt/T sets the resolution limit in uii to distinguish 
two neighboring frequency clusters. A group of oscilla- 
tors is determined to be members of a synchronized clus- 
ter only if all the members (i = n, n + 1, . . . , n + m + 1) 
have the same value of u>i as the neighboring sites within 
some tolerance \uji — u>i\ < i] (£ = i — 1 and i + 1). Since 
T^ 1 is O (l0~ 4 ) the tolerance rj is set to 10~ 3 . 



C. Comparison of RG and simulations. 

While we are primarily interested in correctly captur- 
ing the statistical characteristics through RG, the local 
validity of our procedure can be established by compar- 
ing the frequency cluster structure that it predicts with 
results from the simulations. Figure [3] plots these com- 
parisons over a representative sample of 150 oscillators 
from chains ofN= 10000 at 3 different coupling distribu- 
tion widths. For all 3 values of n, there are clusters with 
large frequency oscillators interspersed. As expected, the 
characteristic size of these clusters grows with fi. There 
is excellent agreement of the cluster sizes, locations, and 
frequencies between the two approaches. This provides 
some confidence that the essential mechanisms are cap- 
tured by the RG procedure. 

Analysis of these and similar plots, shows that the 
agreement of the cluster frequencies becomes more ac- 
curate with increasing cluster size. The reason for this 
result is the smaller boundary-to-bulk ratio of larger clus- 
ters, which is more apparent if Eqs. ([5]) is rewritten as 

9i=c>i + sin (0 t -_i - 9t) + — sin {9 i+1 - 6,) (28) 

rrii m i 



As clusters become larger, the importance of the coupling 
to the rest of the chain is weighted down by a 1/m fac- 
tor. In the fast oscillator decimation, this is manifested 
by the appearance of this factor in the frequency correc- 
tion. Also, when a cluster is built up by a strong coupling 
decimation step, any frequency renormalizations of the 
component oscillators by prior fast oscillation decimation 
steps cancel when the weighted sum of the frequencies is 
calculated to give the frequency of the new cluster. Thus 
the cluster frequencies u)i are given to a good approxima- 
tion by the average of the constituent bare u> values, plus 
the small correction due to the fast oscillator decimation 
step that yields the final cluster. 

Investigating the frequency distribution of clusters of 
a given size provides further insight into the validity of 
this RG procedure. This is shown in Fig. 21 where for 
each cluster size m from 1 to 4, the RG predictions are 
compared with the results from simulations using three 
values of the parameter fx. The two methods are in ex- 
cellent agreement for the larger values of m and also for 
all m at small values of /i. 

The comparisons in Fig. 0] show a small but system- 
atic underestimation of the number of smallest sized clus- 
ters (m = 1 and m = 2) near zero frequency and large 
coupling [panels (b) and (c) and less so in in panels 
(e) and (f)]. The rms difference between the RG gen- 
erated and numerically calculated curves ranges between 
0.7— 2.2-10 -4 , with the worst agreement in panel (b), and 
best agreement in panel (1) . The region of the greatest 
discrepancy, small m and w, is in fact where the fast os- 
cillator decimation is expected to be least accurate, since 
by the time the RG sweep reaches these frequencies, the 
distributions of remaining frequencies and couplings are 
no longer expected to be wide so that the approximations 
based on strong randomness are not as accurate. How- 
ever the disagreement decreases with increasing m, be- 
cause, as described above, the frequencies of larger clus- 
ters are insensitive to the corrections from neighbors. 



IV. THE UNSYNCHRONIZED PHASE OF AN 
OSCILLATOR CHAIN 

Let us now use the results of the RG and the simula- 
tions to discuss the statistical properties of the frequency 
clustering in the one dimensional oscillator chain with 
strong randomness. 



A. Cluster size distribution 

The main measurable physical quantity that arises in 
our analysis is the distribution of cluster sizes. In Fig. [5] 
we show the distribution of cluster sizes for several widths 
\x of the coupling distribution from the RG and simula- 
tions. The RG was performed on one single realization of 
10 6 oscillators, while the simulations were performed on 
100 realizations of 10 4 oscillator chains. Since the char- 
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FIG. 3: Comparisons between cluster structure of dynamic frequencies found in numerical simulations of a 10000-oscillator 
chain with those predicted by RG over a representative number of lattice sites. The top row (a-c) plots the dynamic frequencies 
found for 500 oscillators in the simulations and demonstrate how some oscillators maintain large frequencies. The dashed boxes 
in the centers of (a-c) encircle 50 oscillators magnified in panels (d-f ) , respectively. In (d-f ) average frequencies from simulations 
given by Eq. (|27p are denoted by solid squares, while RG predictions are shown as open circles. The differences between RG 
and simulations are shown in panels (g-i). Each column shows results for a different coupling width: /i = 1.25 (a,d,g), /i = 3.75 
(b,e,h), and \x = 7.5 (c,f,i). 



acteristic cluster size is much smaller then 10 4 , we com- 
bine results from all 100 realizations and analyze them 
as a single chain of size 10 6 . As can be seen in Fig. [5J 
good agreement exists between the RG and the simula- 
tions. The apparent discrepancies at the largest cluster 
sizes are due to statistical fluctuations resulting from the 
small number of such clusters present in the 10 6 oscillator 
ensemble. 

For a given value of the coupling width parameter //, 
the probability of finding a cluster is expected to fall off 
rapidly with increasing cluster size. The linear scaling 
in the semi-log plot suggests that distribution of cluster 
sizes has the form of an exponential: P(n) oc exp (— n/£). 
In Fig. [5] we exhibit the linear fits to the RG data. The 
fits were made between twice the characteristic size £ and 
the cluster size at which there are just 20 occurrences of 
clusters of that given size, since for large clusters the 
scatter becomes too large. 



B. Dependence of £ on disorder strength 

The function is a comprehensive quantity that 

characterizes the physics of the chain in the unsynchro- 
nized regime, because [i is the only control parameter in 
the chain, and £ is the only quantity that characterizes 



the cluster statistics since the distributions of cluster sizes 
are exponential. We plot the function in Fig. [6j We 
find good agreement between the RG procedure (open 
circles) and simulation (solid squares) until /i«5, after 
which the two results begin to deviate from each other. 
The small, but growing discrepancy beyond /i = 5 sug- 
gests there is an additional mechanism in the physics of 
cluster formation that becomes important at larger and 
is not taken into account by the decimation steps of the 
RG procedure. 

A physical effect that may not be adequately included 
in our RG procedure is the accumulation of phase strain 
as a cluster is built up through successive strong coupling 
decimation steps. We can understand this effect based on 
the arguments of Ref. [H[ . If we neglect the interactions 
of a synchronized cluster with the rest of the chain, its 
frequency is the average frequency (uj) of the constituent 
oscillators. Furthermore, forming the sum of Eqs. © 
from one end of the cluster (here labeled 1) to the nth 
oscillator in the cluster gives 

n 

X n = - (u)) = K n sin(0„ +1 - 9 n ), (29) 

i=i 

since the internal interaction terms cancel. We refer to 
this condition for the cluster's existence as a strain check. 




FIG. 4: Number of clusters of a particular frequency divided by the total length of the chain (TV = 10 6 ). The plot compares the 
RG prediction (solid lines) and the simulation results (open diamonds) for 3 coupling widths and 4 cluster sizes. Each panel 
has a fixed coupling width and cluster size. Frequencies, on the hortizontal axes, are grouped into bins of size Aui = 0.1. The 
rows are each a different cluster size; m — 1 (a-c), 2 (d-f), 3 (g-i), and 4 (j-1) from top to bottom. Each column is a different 
coupling width: /i = 1.25 left, /i = 3.75 center, and /i = 7.5 right. 



The condition reads: \X n \ < K n for all n in the clus- 
ter. Note that X n is the sum of random numbers, so 
that as the size of the cluster increases, the range of val- 
ues of \X n \ typically grows, and increasing phase strains 
\0 n+ i — 6 n \ are needed to keep the various parts of the 
cluster synchronized. If the strain becomes so large that 
some \9 n +i —0 n \ become tt/2, the cluster will break. Also 
note that if two clusters are combined through a strong 
coupling decimation step to form a putative larger clus- 
ter, the mean frequency (ui) to be used in the sum of 
Eq. changes, and the condition \X n \ < K n may be 
violated for a bond in the interior of the two subclusters. 
This means that the combined set of oscillators are not in 
fact synchronized, and should be broken into smaller syn- 



chronized clusters in the RG procedure. We proceed by 
making the assumption that the oscillators form just two 
synchronized clusters, and identify the break between the 
two clusters as where the condition \X n \ < K n is violated 
in the presumed cluster. Thus if X n exceeds K n at some 
n, we do not allow the cluster to form. Instead, we re- 
form the bare oscillators into two different sub-clusters, 
such that the first sub-cluster is formed out of the bare 
oscillators 1 through n and the second is formed form 
the remaining bare oscillators. The RG then proceeds as 
normal. The new sub-clusters may be subjected to either 
of the two decimation steps later in the sweep, or in the 
next sweep if more are required. 

The obtained with the version of the RG which in- 
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FIG. 5: Number of clusters over a range of sizes. The plot 
compares RG predictions (open circles) with numerical inte- 
grations of Eqs. (f5]) (solid squares). The curves are spread 
out for clarity by multiplying each successive data set as /j, 
increases by 0.1 (a linear shift in the semilog plot). The solid 
lines denote the fits to the RG data over the domain described 
in the text. 



eluded the strain check are shown by diamonds in Fig. [6] 
The change in the results takes place in the right direction 
(decreasing £), but with a systematic over-estimation. 
All three sets of data, when plotted on a log-log scale fol- 
low approximately straight lines for [i > 0.5, so £(//) ap- 
pears to become a power law in this range. The straight- 
line fits to the data were made for > 0.625. The dashed 
curve in Fig. [5] is the fit to the simulation results, sug- 
gesting the behavior £ oc /i - 46 while the fits to the RG 
data with and without the strain check give £ oc fi 0AS 
and /i 48 . 



C. Estimate of £(/Lt) 

Given the simplicity of the renormalization group ap- 
proach, one may wonder whether we can capture the 
behavior of the oscillator chain with a simple one-loop 
argument, based on the role of weak couplings in the 
limitation of cluster sizes. 

Indeed, the small clusters are likely to form between 
weak couplings, and since the distribution of separations 
between such weak couplings follows Poisson statistics, 
an exponential distribution arises naturally. Using this 
hypothesis the probability of finding a cluster of size n 
(i.e., a cluster delimited by two weak couplings n units 
apart) is P(n, p) = C(l - p) n = Ce~ n ^ where £ = 
— l/ln(l — p) and p is a probability that any randomly 
chosen coupling is weak. We make an estimate of this p 
based on bare distributions of u> and K. Similarly to the 
RG, we define a coupling to be weak if 2K < \oji — u) r \. 
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FIG. 6: Characteristic cluster length £ given by the fits to 
data in Fig. [5] (a): Results of simulation (solid squares), 
the numerical RG of Sec. [TT] (open circles) and an enhanced 
numerical RG with strain check (diamonds) described in 
Sec. IIVBI The solid line is a simple analytic prediction 
£(/i) = — l/log(l — p(/i)) where p{fi) is defined in Eq. (I31[) . 
The dashed line is a power law fit to the solid squares for 
/i > 0.625 and has an exponent of 0.46; similar fits to the dia- 
monds and circles give exponents of 0.48 and 0.48 respectively, 
(b): Log- log plot of the same data excluding the analytic pre- 
diction. 



For each such bond, the probability that the required 
ratio is > 1 is given by 

r°a / r°° \ ( r s / 2 \ 

p = 2 I / g(cj)g(u + S) duj J x / G{K) dK dS 

(30) 

The integral of g(u>)g(ui + 5) is the probability that for 
each of the two oscillators connected by this bond, u>i — 
io r = S. For each such frequency difference, all the bonds 
with K < | S |/2 are considered weak. Finally, all possible 
values of 5 are integrated in the outer integral. After 
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performing the inner integrals over ui and K we obtain 
i /"oo / 5 \ 

tM = *s] ITP arctan {^) dS - (31) 

This result is compared with both the RG predictions 
and numerical solutions in Fig. [6l The expression gives a 
good description of the data for /j->0, but overestimates 
the characteristic cluster size for larger values of /i. The 
renormalization group method, which does not treat all 
scales of frequencies and couplings at once, clearly pro- 
vides better predictions then an estimation based on bare 
distributions. 



V. CONCLUSIONS 

In this paper we investigated the collective behavior of 
a one-dimensional chain of coupled non-linear oscillators 
with random frequencies and nearest-neighbor couplings. 
For this purpose we developed a real space renormaliza- 
tion group approach, which is expected to be reliable in 
analyzing chains with large disorder. As we report above, 
our RG approach did an excellent job in capturing both 
the dynamics of individual oscillators, and the large scale 
behavior of the chain. 

The disordered oscillator chain has the potential of es- 
tablishing macroscopic order, which is exhibited by a fi- 
nite fraction of all oscillators in the system moving in 
accord with the same frequency. Indeed in systems with 
higher connectivities, such global synchronization may 
occur, but in the limit of short ranged interactions, the 
macroscopic order is stymied by fluctuations. This dy- 
namical behavior is reminiscent of critical phenomena in 
equilibrium statistical mechanics, where for each system 
there exists a lower critical dimension, and only dimen- 
sionality higher than that allows macroscopic ordering. 
In our analysis, we concentrated on the quantities that 
most directly capture the collective aspects of the chain's 
behavior, and therefore also the competition between the 
interactions which seek to establish a synchronized phase, 
and the disorder-induced fluctuations which destroy it. 
This physics is fully contained in the cluster-size distri- 
bution, P(n) ~ e~ n /t W , which we find as a function of 
the interaction-strength tuning parameter, fi. 

The behavior of the average cluster size, acts 
as an effective correlation length for the oscillator chain, 
where we once more think of our system in terms of the 
theory of critical phenomena. Since our oscillator chain 
is below its lower critical dimension, it's phase diagram 
resembles that of the one-dimensional Ising model, with 
\x serving the role of inverse temperature 1/T. The chain 
has an ordered phase, where all oscillators are synchro- 
nized, which is reached in the limit /i — » oo (analogous 
to the Ferromagnetic phase at T — ► 0). This phase, 
however, is described by an unstable fixed point. If an 
RG flow can be extended to the entire parameter space, 
the full-synchronization fixed point flows under RG to 
the stable fixed point found at the non-interacting point, 
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FIG. 7: Fixed points and flow diagram of the oscillator chain. 
The short-range oscillator chain has only two fixed points: the 
unstable fully synchronized fixed point at infinite interactions 
(or zero disorder) — » oo, and the stable unsynchronized 
fixed point. The cross over flow between the two points is 
associated with a correlation length, captured by the average 
cluster size, £(/u). 

fj, = 0, which describes the unsynchronized phase (anal- 
ogous to the paramagnetic fixed point at T — > oo). The 
putative flow diagram of the chain is shown in Fig. [7] 
The flow from the fully-synchronized point to the un- 
synchronized fi = point is indeed associated with a 
length scale, whose meaning is the coarse-graining length 
at which the chain seems fully unsynchronized, i.e., con- 
sisting of single free oscillators. This length is the size of 
synchronized clusters in the original chain, £(/z). 



Since £ (fi) serves as an effective correlation length, it may 
diverge in a universal fashion in the limit /i — > oo. While 
our RG method can only access moderate values of fi, we 
cautiously argue that we can access the critical region, 
as can be seen in Fig. [5l From our results we obtain the 
critical exponent v of the synchronization transition: 

£(/z) ~\i v v ~ 0.48 (32) 

Note that this behavior is not analogous to the divergence 
of the Ising chain correlation length £ ~ In y . 

The method by which we obtained the above results 
constitutes the main achievement of this work. We de- 
veloped a real space renormalization group scheme that 
successfully predicts the detailed behavior of a nonlinear 
oscillator chain with short-range interactions. We have 
implemented the RG scheme numerically on a chain of 
10 6 oscillators and compared its results with essentially 
exact simulations of the model equations, Eqs. ([3]) (see 
Sec. IIIip . The agreement is very good, and for some char- 
acteristics, such as the cluster frequency distributions at 
low frequencies, it improves with larger cluster sizes. This 
indicates particularly that our local RG scheme success- 
fully captures the non-trivial collective features of the 
chain expressed in large synchronized clusters. 

One practical implication of our analysis may be found 
when considering coupled laser arrays, where only limited 
connectivity can be achieved. Diode lasers, for example, 
are manufactured in one-dimensional bars that allow for 
evanescent coupling through overlapping electric fields. 
Due to the exponential fall-off in the electric field enve- 
lope with distance, the coupling between well-behaved 
diodes can only be short ranged. The tools developed 
here could be used to study the extent of coherence that 
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such systems exhibit. A potentially promising area of 
current research is constructing grids of diodes coupled 
in 2-dimensions by stacking 1-dimensional bars. Our RG 
technique could possibly also be extended to such a sys- 
tem, and help clarify both its properties, and whether it 
can establish macroscopic coherence. 

Our work suggests several other directions for future 
research. First, we may ask how can one improve the 
applicability of our RG scheme to larger values of the in- 
teraction parameter \i. This may be pursued by focusing 
on the concept of intra-cluster strain, which is roughly 
the variance of phase differences between neighboring os- 
cillators within a cluster. As interaction increases, cluster 
sizes grow, and bigger strain is needed to accommodate 
the spread of frequencies within a cluster (cf. Sec. IIV Bl) . 
In our current scheme we made the first attempt to take 
this into account with the strain check based on the argu- 
ment of [HI], a condition expressed by Eq. (|29|) . This led 
to some improvement - for example, was modified 
in the right direction, but this modification was over- 
estimated. The strain check may be augmented in the 
future. 

A second natural question is the degree of universal- 
ity in the oscillator chain. Our result for the correlation 
length critical exponent v, Eq. (|32|) . was obtained using 



Lorentzian distributions. It is possible that the exponent 
depends on the nature of the tails of the distributions 
used; Lorentzians, for instance, do not have a variance. 
We intend to investigate the universality of v alongside 
investigating the applicability of our RG scheme to dis- 
tributions with narrower tails. 

A more general question concerns the development of 
an RG scheme and the analysis of synchronization in 
short-range networks at higher dimensions, and espe- 
cially at two dimensions, which is very close to the lower 
critical dimension: Refs. [13, HI] suggest that in two di- 
mensions there may be a transition to a macroscopically 
synchronized phase, or that d = 2 is the lower critical 
dimension giving rise to interesting anomalous behavior 
(see also [13]). The development of a decimation proce- 
dure at dimensions higher than one is quite challenging 
because of the higher connectivity of the system, and the 
possibility that real-space local decimation steps change 
the geometry of the system. 
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ous discussions and Tony Lee for carefully checking the 
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National Science Foundation under Grant No. DMR- 
0314069 for funding this work. 
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